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I from AGS to RHIC decreases the value of the chemical potential in the Hadron- 
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I QGP phase diagram. Thus, the formation of QGP at RHIC energies in central 

^ , collisions may be expected to occur at very small values of the chemical poten- 

I tial. Using many different relativistic mean-field hadronic models (RMF) at this 

.!_H ' regime we show that the critical temperature for the Hadron-QGP transition is 

X: 

H hadronic model independent. We have traced back the reason for this and con- 

elude that it comes from the fact that the QGP entropy is much larger than the 
hadronic entropy obtained in all the RMF models. We also find that almost all of 
these models present a strong entropy enhancement in the hadronic sector com- 
ing from the baryonic phase transition to a nucleon-antinucleon plasma. This 
result is in agreement with the recent data obtained in the STAR collaboration 
at RHIC where it was found a rich proton-antiproton matter. 
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I. INTRODUCTION 



The understanding of nuclear matter under extreme conditions is a crucial and indispensable 
aim of nuclear and stellar physics, especially with much more experimental information to come 
with RHIC and Alice/LHC accelerators. It is now widely believed that at sufficiently high 
energy, in heavy- ion collisions, a central hot region is formed This region is commonly 

associated with the existence of a Quark-Gluon Plasma (QGP). Following many prescriptions, 
this hot region expands, cools down and freezes out into hadrons, exhibiting a hadron-gas 
phase transition. In order to observe this transition to the deconfined phase, some observable 
signatures have been proposed as, for example, electromagnetic radiation strangeness 
enhancement 0], and i /ip suppression High-energy heavy- ion collisions data has brought 
still more excitation about the predictions of QCD on this subject. In fact, based on the results 
from different experiments done at CERN-SPS, it is believed that a new state of matter has 
been formed, the quark matter ||T0 |. 

Recently, phase transition signatures at resonance-rich matter in heavy ion collisions at 
RHIC energies has been analyzed in central Au + Au collisions at \/s= 200 A GeV [|n|)IIl]- 
These studies, using a microscopic transport model (UrQMD) calculation, showed the equili- 
bration of hot matter through a. Txfi phase diagram. Confronting these results with RHIC, 
SPS and AGS data, it was shown that for central collisions the rising of the incident energy 
implies a strong evidence for hadronic-QGP phase transition at very small baryon chemical 
potentials and temperatures lying between 150 to 200 MeV. 



Some time ago, using the resonance model, Heinz et al. [T^ , have studied the Hadron-QGP 
phase transition at this regime. They found a strong enhancement of the entropy from the 
hadron to QGP phase. When they included all the resonances and mesons, the entropy of the 
hadron phase increased but a considerable entropy gap still remained between the two phases. 
In this paper we will study this transition in relativistic mean field (RMF) models, where we 
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have not only a repulsive vector interaction but also an attractive scalar interaction, that was 
not considered in the resonance model study . This scalar interaction, as we will show in 
this paper, is quite important in order to study this hadron-QGP transition at low net baryon 
density. At this regime there is no vector field and the pressure is composed only from the 
nucleon thermal gas and from the scalar field pressures. The abrupt fall in the nucleon mass at 
high temperature, seen in many of the RMF hadronic models, is due to the abrupt increasing of 
the scalar field. This effect resembles a phase transition when the system becomes a dilute gas 
of baryons in a sea of baryon-antibaryon pairs, and it is caracterized by a large enhancement of 
the entropy |]n|,|TB|,|TB[ . Thus, our study in RMF models is quite important because the entropy 
enhancement depends strongly on the scalar interaction presented in these type of models. It 
is important to stress, that in all these models the entropy enhancement is maximized at zero 
chemical potential because when we increase the baryon density, the fall of the nucleon mass 
with the temperature and also the increase of the entropy become smaller [|16|. The importance 
of the entropy production in relativistic heavy ions collisions has also been addressed to analyze 
the QGP phase transition, where the entropy density jump at the phase boundary for low net 



baryon density has been pointed out |17 



At the hadronic level, RMF models are able to describe well a number of nuclear phenomena 
through different conceptions of meson-nucleon and meson-meson couplings. Most of them 
follows the basic Walecka model [0] in which improvements have been done in order to better 
fit the finite nuclei data |T^J2yj2T|J2^ . It is known that hadronic models may present phase 



transitions at higher temperatures. For example, in reference [14|, hot nuclear matter was 



extensively studied using the Walecka model regarding the finite temperature behavior for zero 
baryon density (nucleon-antinucleon plasma). The authors explored this model for different 
scalar- vector coupling constants to conclude that, depending on the values taken by those 
constants, the Walecka model may or may not generate a phase transition into a nucleon- 
antinucleon plasma. The order (first or second) of this phase transition depends crucially on 



the scalar coupling constant of the model. 

The question we pose in this paper is whether the phase transition for the hadronic sector 
itself, described from hadronic models, may or not have importance to the QGP phase transi- 
tion. To answer this question we bring to this discussion a representative class of recent RMF 
hadronic models . We have also considered derivative coupling models pO|j2T[| and different 



parametrizations of the non- linear Walecka model such as NLl, NL2, NL3 and NLSH models 
very well known in the literature [^. Most of these models were successfully employed to cal- 
culate nuclear matter bulk properties as well as the spectra of finite nuclei. The aim of our work 
is to present a systematic comparative study of a set of hadronic models, at extreme temper- 
atures and very low net baryon density, and the hadronic-QGP phase transition, emphasizing 
the entropy density behaviour in both phases. The QGP phase is represented by a perturbative 
QCD derivation We studied the hadronic-QGP phase transition for a broad range of 

bag constant B and QCD running coupling constant ag values. The particular case of free 
quarks [og = 0) is also explicited. This kind of phase transition has been investigated by many 
authors p4| , p5| . In particular, aspects of the role played by the nucleon exclude volume have 
also been studied p6[| . 

As a result, we find that practically all the analyzed hadronic models signalize approximately 
the same critical QGP temperature. This finding is quite independent whether the hadronic 
models have or not hadronic phase transition themselves at high temperature. This suggests 
strongly a RMF hadronic model independence for the QGP phase transition in this regime. To 
the best of our knowledge, this is a remarkable and unknown feature of the hadronic models. 
We have traced back the reason for this and conclude that it comes from the fact that the QGP 
entropy (the slope of the pressure versus temperature) is much larger than the hadronic entropy 
obtained in all the hadronic models. Therefore, insofar the entropy through sophisticated 
calorimeters in the future may be a measured quantity, we believe that at very low net baryon 
density, the enhancement of the entropy may become a signature of a QGP. Previous works 



have already addressed the role played by the entropy at the QGP formation |27|j28|JT7|l . We 
have observed a large jump of the hadronic entropy density, seen in almost all of the RMF 
models. An enhancement of the entropy per baryon with the rising of the incident energy was 



also found in recent microscopic transport model (UrQMD) calculations [TT||T^. However, this 
model did not include a QGP phase but can be taken as leading credence for our hadronic 
model phase entropy enhancement signals. This enhancement may favour the hadron-QGP 
phase transition once it decreases the latent heat for the QGP formation. Thus, we expect in 
central collisions with a low net baryon density, the formation of a very rich baryon-antibaryon 
matter in the region of the QGP phase. This results is in agreement with recent data obtained 
in the STAR collaboration at RHIC where in fact it was found a rich proton-antiproton matter 

The outline of this paper is as follows. In Sec. II we present the phase transition calculation. 
In section III we show our results and discussion, followed by the conclusions in section IV. 



II. THE PHASE TRANSITION CALCULATION AT HIGH T AND ^ = 

At low temperatures and baryon densities (around the nuclear matter ground state), renor- 
malization group arguments |3y] show that the QCD running coupling constant is greater than 
one, indicating confinement of quarks and gluons inside hadrons which are the appropriate 
degrees of freedom of the nuclear matter in this regime. In this region, QHD based on effective 
Lagrangians works quite well. The phase transition is supposed to occur at sufficiently high 
temperatures and densities, where in this case the QCD running coupling constant becomes 
smaller than one, suggesting deconfinement of quarks and gluons. If one assumes the validity of 
QHD far beyond the normal nuclear matter ground state, the transition can be characterized by 
the process of hadrons loosing their identities, and quarks and gluons becoming the elementary 
degrees of freedom. 
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In order to see the hadronic-quark-gluon phase transition we need distinct models for the 
two different phases of the baryonic matter. We describe the hadronic phase (H) by the models 
of Ref. Ill9| , p0| , ^ , p2| in the Hartree approximation extended to finite temperature at zero net 
baryon density. For the quark-gluon plasma phase, the equation of state is given given by 
||T3| , |23[] . A study of the hadron-QGP phase diagram for finite baryonic density of the models 
presented in Ref. [|I9[ has been done in Ref. |]3T|. The Gibbs criteria used in the analysis of the 



phase transition for zero density is given by uniTc) = nqcpiTc) = and Ph(Tc) = pqcpiTc) . 

The functional, used for the hadronic phase, originates from an effective Lagrangian, with 
additional physical constraints ensuring QCD symmetries |T3. The pressure is 



gl \2 3! M 4\ J 



2gl \ 2 M2 3! M3 



where $ and W are the scalar and vector meson fields respectively , E*{k) = (k^ + M*^)^/^ and 
f±{k,T) stand for the Fermi-Dirac distribution for baryons (antibaryons). The baryon density 
Pb is defined by pB = J cPk [/+(/c,T) — /_(A;,T)]. Thus, in the zero net baryon density 
regime we have the same number of baryons and anti-baryons. 

The constants (see Table I) are already adjusted to reproduce some of the bulk properties 
of the nuclear matter ground state [|l^]. Here Wl stands for the usual Walecka model. CI and 
Ql are nonlinear a models, having scalar field self-interaction added to the Wl model in order 
to improve some nuclear matter bulk properties of Wl. Most of the developed hadronic models 
has the CI structure, given by cubic and quartic scalar field self-interaction . In this aspect 
it is a very representative model category. The scalar cubic and quartic self-interaction are 
usually claimed to simulate three-body and four-body forces effects respectively. The others 
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as Q2, Gl and G2 are models representing further improvements by including different scalar- 
scalar, vector-vector and scalar-vector fields couplings. We have also considered for the hadronic 
phase derivative coupling models |2^j2l| and different parametrizations of the nonlinear Boguta- 



Bodmer model gg, NLl, NL2, NL3 and NLSH models fSg. 

For the quark-gluon plasma phase, at high temperature, the pressure and the density is 
given by the expansion 



PQGp{^^q,Tq) = — — 1 



45 V 47r 
47r2 







60 











^ 50a, 



2l7r 

(2) 



PQ- = ^iv,Lr,^ + ^)(i-^), (3) 



where B is the bag constant, Nf is the number of flavors, a<j is the QCD running coupling con- 
stant, depending on the quark-gluon plasma temperature Tq and the quark chemical potential 
/ig through the first order perturbative expression 

2A^.^ -^-^ 



a. 



4vr{(ll - ^) ln[(0.8/.; + 15.622T,^) / K^]] . (4) 



At the zero density regime we are interested we set /i^ = 0. 

The model for the QGP phase which we use in this work has some parameters such as A, 
B and Nj. As a first approximation, we examine the simplest case with Nf = 2 (quarks u and 
d only). The QCD scale parameter A is fixed at 200 MeV, consistent with the current data 
set. We will use two different bag constants in our analysis: B^^^ = 174 MeV and B^^^ = 238 
MeV, corresponding to 5 = 119 MeV fm~'^ and B = 418 MeV fm~^ respectively. These are 
the limiting values of the broad range of values used in the literature. 
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III. RESULTS AND DISCUSSION 



Figure 1 shows the effective nucieon mass ratio m* vs. T for the studied models. We see 
that ZM, ZM3, CI and G2 models do not present signals of a first order phase transition in 
opposite to the others. Then, not all of the hadronic models considered here (twelve in total) 
have the hot nucieon- ant inucleon phase transition in the hadronic sector, manifested by the 
large decrease of the effective nucieon mass M*. 

By using the Gibbs criteria of phase equilibrium Pi = P2i l^i = 1^2 — ^ and Ti — T2 we have 
obtained the critical temperatures for NLl ( 188MeV ), NL3 ( 191MeV ), NLSH ( 192MeV 
),W1 ( 186MeV ), gi ( 192MeV ), Q2 ( 192MeV) and Gl ( 191MeV ) models. The other 
models studied in this paper do not present first order hadronic phase transition. The critical 
temperatures for this hadron nucleon-antinucleon phase transition for hot nuclear matter at the 
hadronic sector itself are remarkably close to the critical temperature of the hadron-QGP phase 
transition. This hadron nucleon-antinucleon phase transition where the effective nucieon mass 
becomes very small can be seen as the Quantum Hadrodynamics version of the chiral phase 
transition. It is interesting that in recent lattice calculations the chiral and the deconfinement 
transitions coincide and in our calculation both transitons seem also to appear at the same 
temperature depending on the value of the bag constant B we choose. 

We present in figure 2 the pressure for hadronic-QGP phases (in units of the Stefan- 
Boltzmann pressure) as a function of T. In this figure we have used the extreme values 
B^^* = 174 MeV and B^^* = 238 MeV. It is to remark how the QGP critical temperatures are 
close for all the hadronic models. Prom this figure we see the critical temperatures very close 
to 151 MeV for = 174 MeV and around 198 MeV for = 238 MeV . The reason lies in 
the slope for the QGP pressure. This slope is nothing but the entropy density, that we show in 
figure 3. Notice how far, by a factor around five ( for cts = ) and around three (for as ^ ), 
the entropy amounts are different in both phases. It causes the hadronic model independence 
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for the hadronic-QGP critical temperature. 

When the hadronic model itself presents the nucleon-antinucleon phase transition it in- 
creases substantially the entropy density as can be seen in figure 3. However, this enhancement 
is not enough when compared with the QGP entropy. This shows a discontinuity in the en- 
tropy signalizing a first order hadron-QGP phase transition, once we are committed with only 
nucleon-antinucleon in the hadronic phase. Let us also note that the QGP entropy density is 
independent of the bag constant B since it is given by {Sqgp + Pqgp )/T. Note how, in figure 
2 for different bag constants, the two slopes are the same for pqgp ■ However, as pointed out 
above, the bag constant has effect on the hadronic-QGP value of Tc . Therefore, if one thinks 
that the enhancement of the hadronic entropy density would help the QGP formation, only 
values of S^/^ greater than 200MeV should be considered. 

We have then included, in an ad hoc fashion, a thermal pion gas contribution to the hadronic 
phase. By analyzing the behaviour of S/p as a function of T with and without pions we conclude 
that the change in the critical phase transition temperature is very small. When we have only 
nucleons and pions in the hadronic phase, the hadronic-QGP phase transition at very low net 
baryon density takes place for approximately the same value of Tc in all the studied hadronic 
models. The hadronic-QGP phase transition is essentially given by the QGP phase and lies 
very close to the temperature where pqgp crosses zero from below. 

The chiral transition version we found in the hadronic models can be identified by a strong 
decreasing in the effective baryon mass and an abrupt increasing of the entropy density in the 
hadronic phase ( see Fig. 1 and Fig. 3) . This large hadronic jump of the entropy density favours 
the hadron-QGP phase transition because it requires a smaller latent heat at the transition 
compared with the models that do not show that pure hadronic nucleon-antinucleon plasma 
phase. Thus, from our results, we would expect at very low net baryon density a formation of 
a very rich baryon-antibaryon matter just before the QGP phase. 

It is important to stress that, even with the formation of this rich baryon plasma, the 
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entropy enhancement is not enough yet to be compared with the QGP entropy density and 
still requires a large latent heat at the phase transition. The latent heat and entropy at the 
transition, as it has been discussed before on the /i = central rapidity region p3| , ^^ 



may be seen by a large change in the specific volume that slows down the time scale for the 
conversion of QGP into the hadron matter: a large entropy discontinuity implies a long lifetime 
of the plasma. This point itself, the entropy discontinuity, has deserved intense work and 
analysis from many authors In fact it addresses direct the order of the hadronic-QGP 

phase transition. We are aware of the limitations of our analysis. The results are dependent 
on the QGP equation of state we are using, obtained by a perturbative expansion in which 
we treat the quarks as massless. Note in figure 3 how the entropy density decreases with the 
increase of the running coupling constant that at the phase transition is around 0.4 to 0.6 . 
This suggests that a non perturbative treatment for the QGP equation of state would reduce 
the entropy density. Indeed, recent lattice QCD results for the pressure and entropy density 
36| , |37[| have also obtained a strong entropy enhancement close to the phase transition exactly 



at the order of our critical hadronic entropy density (Figure 3). These lattice results were also 
reproduced by a non perturbative purely gluonic QCD calculation [^]. Therefore, we should 
say that the order of the hadron-QGP phase transition is still an open point in our work. 



IV. CONCLUSIONS 

Using many different RMF models for the hadronic phase and a perturbative QCD equation 
of state for the QGP, we study the hadron-QGP phase transition at zero net baryon density 
and high temperature. We show that the critical temperature for this transition is hadronic 
model independent. We have traced back the reason for this and conclude that it comes from 
that the QGP entropy (the slope of the pressure versus T) is much larger than the hadronic 
entropy obtained in all the RMF models. This finding is quite independent whether the hadronic 
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models have or not a hadronic phase transition at high temperature when the system becomes 
a dilute gas of baryons in a sea of baryon-antibaryon pairs. 

The hadronic-QGP phase transition at zero net baryonic density takes place for approxi- 
mately the same value of Tc in all the studied hadronic models. Among the studied hadronic 
models most of them showed a pure hadronic phase transition, the Quantum Hadrodynamics 
version of the chiral phase transition, seen by a strong decreasing in the effective baryon masses 
and an abrupt increasing of the entropy density in the hadronic phase. This fact by itself 
does not change the hadronic-QGP value of Tc . However, it is important to remark that the 
latent heat given by Tc {Sqgp — Sh) connecting both phases is much lower (but still not small, 
Sqgp/Sh = 3 with a non zero ag, see figure 3) for the hadronic models which present this 
hadronic phase transition to a baryon-anti-baryon plasma, favouring the Hadron-QGP tran- 
sition. This result enlarges the knowledge of this transition regarding RMF models and may 
reinforce different model calculations conjecturing the enhancement of the entropy as a signa- 
ture for a hadronic-QGP phase transition at low net baryonic density. This signature it is quite 
important nowadays with the underway experiments, such as the collaboration STAR at RHIC, 
because recent theoretical studies showed that for central collisions the rising of the incident 
energy from AGS to RHIC decreases the value of the chemical potential in the Hadron-QGP 
phase diagram. Thus, the formation of QGP at RHIC energies may be expected to occur at 
very small values of the chemical potential, exactly the regime we are analyzing in this pa- 
per, were we expect a jump on the hadronic entropy density due to the formation of a rich 
baryon-antibaryon matter. Finally, we stress the importance of the direct measure through 
calorimeters of this quantity. 
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TABLES 



TABLE I. Model parameters, taken from Ref. [3]. 



Model 


Wl 


CI 


Ql 


Q2 


Gl 


G2 


rUs/M 


0.60305 


0.53874 


0.53735 


0.54268 


0.53963 


0.55410 




0.93797 


0.77756 


0.81024 


0.78661 


0.78532 


0.83522 




1.13652 


0.98486 


1.02125 


0.97202 


0.96512 


1.01560 


Vi 




0.29577 






0.07060 


0.64992 


m 










-0.96161 


0.10975 






1.6698 


1.6582 


1.7424 


2.2067 


3.2467 


K4 






-6.6045 


-8.4836 


-10.090 


0.63152 


Co 








-1.7750 


3.5249 


2.6416 
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FIGURES 
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FIG. 1. Effective nucleon mass as a function of the temperature for all the hadronic models. 




FIG. 2. The pressure for hadronic-QGP phases (in units of the Stefan-Boltzmann pressure) as 
function of T. 
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FIG. 3. Entropy density as a function of T for all the hadronic models and for the QGP plasma 
the QGP entropy density is divided by 3 ) . 
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